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ABSTRACT 

Aims. The questions of how strong magnetic fields can be stored in rotating stellar radiative zones without being subjected to pinch- 
type instabilities and how much radial mixing is produced if the fields are unstable are addressed. 

Methods. Linear equations are derived for weak disturbances of magnetic and velocity fields which are global in horizontal dimen- 
sions but short-scaled in radius. The equations are solved to evaluate the stability of toroidal field patterns with one or two latitudinal 
belts under the influence of a rigid basic rotation. Hydrodynamic stability of latitudinal differential rotation is also considered. 
Results. The magnetic instability is essentially three-dimensional. It does not exist in a 2D formulation with strictly horizontal distur- 
bances on decoupled spherical shells. Only stable (magnetically modified) r-modes are found in this case. The instability recovers in 
3D. The most rapidly growing modes for the Sun have radial scales smaller than 1 Mm. The finite thermal conductivity makes a strong 
destabilizing effect. The marginal field strength for the onset of the instability in the upper part of the solar radiative zone is about 
600 G. The toroidal field can only slightly exceed this critical value for otherwise the radial mixing produced by the instability would 
be too strong to be compatible with the observed lithium abundance. Also the threshold for hydrodynamic instability of differential 
rotation which exists in 2D is lowered in 3D. When radial displacements are included, the value of 28% for critical shear is reduced 
to 21%. 
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1. Introduction 

The old problem of hydromagnetic stability of stellar radia- 
tive zones attains a renewed interest in relation with Ap stars 
magnetism (Braithwaite & Spruit I2004I I. the solar tachocline 
(Gilman 2005 ) and transport processes in stably stratified stel- 
lar interiors (Barnes et al. U999b . 

Which fields the stellar radiative cores can possess is rather 
uncertain. The resistive decay in radiation zones is so slow that 
primordial fields can be stored there (Cowling 1945 ). Whether 
the fields of 10 5 G which can influence g-modes of solar oscilla- 
tions (Rashba et al. 2006 ) or still stronger fields which can cause 
neutrino oscillations (Burgess et al. 120031 1 can indeed survive in- 
side the Sun mainly depends on their stability. 

Among several instabilities to which the fields can be sub- 
jected (Acheson |1978l l, the current-driven (pinch-type) instabil- 
ity of toroidal fields (Tayler 1973) is probably the most relevant 
one because it proceeds via almost horizontal displacements. 
The radial motions in radiative cores are strongly suppressed 
by buoyancy. Watson (1981) estimated the ratio of radial (u r ) 
to horizontal (ug) velocities of slow (subsonic) motions in rotat- 
ing stars as u r /ug ~ Q. 2 /N 2 , where Q is the basic angular ve- 
locity and N the buoyancy (Brunt-Vaisala) frequency. This ratio 
is small in radiative cores (Fig. []]). If the radial velocities are 
completely neglected, the stability analysis can be done in 2D 
approximation with decoupled spherical shells (Watson 1981; 
Gilman & Fox I1997I ). There is, however, some radial motion 
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still excited. How much radial mixing the Tayler (119731 ) insta- 
bility produces is not well-known so far. As the radial mixing is 
relevant to the transport of light elements, a theory of the mixing 
compared with the observed abundances can help to restrict the 
amplitudes of the internal magnetic fields (Barnes et al. 1999 ). 
In the present paper the vertical mixing produced by the Tayler 
instability is estimated and then used to evaluate the upper limit 
on the magnetic field amplitude. 

In this paper, the equations governing linear stability of 
toroidal magnetic fields in differentially rotating radiation zone 
are derived. They are solved for two latitudinal profiles of the 
toroidal field with one and with two belts in latitude but only 
for rigid rotation. The unstable modes are expected to have the 
longest possible horizontal scales (Spruit 1999). Accordingly, 
our equations are global in horizontal dimensions. They are, 
however, local in radius, i.e. the radial scales are assumed short, 
kr » 1 (k is radial wave number). The computations con- 
firm that the most rapidly growing modes have kr ~ 10 3 but 
they are global in latitude. The derived equations reproduce the 
2D approximation as a special limit. It is shown with an ex- 
actly solvable case that for rigid rotation the Tayler instability 
is missing in the 2D approximation. It recovers, however, in 3D 
case. Finite diffusion is found important for the instability. The 
minimum field producing the instability is strongly reduced by 
allowance for finite thermal conductivity. This field amplitude 
is about 600 G for the upper part of the solar radiation zone. 
Considering the light elements transport by the Tayler instabil- 
ity we find that the field strength can be only slightly above the 
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Fig. 1. The buoyancy frequency (fj) in the upper part of the ra- 
diative core of the Sun after the model by Stix & Skaley (1990 ). 
The convection zone of the model includes the overshoot layer 
so that N/Q.Q is large immediately beneath the convection zone. 



marginal value. Otherwise, the mixing would not be compatible 
with the observed lithium abundance. 



2. The model 

2.1. Background state and basic assumptions 

The stability of rotating radiation zone of a star containing mag- 
netic field is considered. The field is assumed axisymmetric and 
purely toroidal. Hence, it can be written as 



B - e$r sin 9 -y^op Qa (r, 9) 



(1) 



in terms of the Alfven angular frequency O a . In this equation, 
p is the density, r, 9, <p are the usual spherical coordinates and 
#0 it the azimuthal unit vector. Equation ([T} automatically en- 
sures that the toroidal field component vanishes - as it must - at 
the rotation axis. Centrifugal and magnetic forces are assumed 
small compared to gravity, g/r » Q. 2 and g/r » Q^. Deviation 
of the fluid stratification from spherical symmetry can thus be 
neglected. 

The stabilizing effect of a subadiabatic stratification of the 
radiative core is characterized by the buoyancy frequency, 



N 2 = -^-, 
C p 8r 



(2) 



where s — C v In (P/p Y ) is the specific entropy of ideal gas. The 
frequency N is very large in the radiative cores of not too rapidly 
rotating stars like the Sun (N » Q., N » Q. A , see Fig.Q]). 

The larger is N the more the radial displacements are op- 
posed by the buoyancy force. Radial velocities should therefore 
be small. They are often neglected in stability analysis what 
might be dangerous as certain instabilities may even be sup- 
pressed by the neglect (we shall see later how this indeed hap- 
pens for the Tayler instability). 

The consequences of the neglect of the radial velocity pertur- 
bations, u' r , can be seen from the expression for (divergence-free) 
velocity, u', in spherical geometry in terms of the two scalar po- 
tentials for the poloidal, P v , and toroidal, P v , flows 

„' = fl£p ggf 1 dTv | ° 2P A | 
r 2 v r \ sin 9 8if> drd9 ) 
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(cf. Chandrasekhar 1 1961ft . Here and in the following distur- 
bances are marked by dashes. For zero radial velocity, the 
whole class of disturbances corresponding to the poloidal flow 
vanishes. Then the remaining toroidal flow can produce only 
toroidal magnetic disturbances so that in the expression for the 
magnetic field perturbations (B'), 



B' = B -L£p ee { 1 dTm i dlPm 
r 2 ™ r \sin0 dcp 8r89 
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the poloidal field potential becomes zero, P m = 0. Reducing the 
class of disturbances may switch-off some instabilities. 

It can be seen from Eq. (01 that the horizontal part of the 
poloidal flow can remain unchanged when the radial velocity is 
reduced and the radial scale of the flow is reduced proportion- 
ally. The disturbances can thus avoid the stabilizing effect of the 
stratification by decreasing their radial scale. 

Our stability analysis is local in the radial dimension, i.e. 
we use Fourier modes exp(iAr) with kr » 1 . It will be confirmed 
that the most unstable modes do indeed prefer short radial scales. 
The analysis remains, however, global in horizontal dimensions. 

Instabilities of toroidal fields or differential rotation proceed 
via not compressive disturbances. Characteristic growth rates of 
the instabilities are small compared to p-modes frequencies. In 
the short-wave approximation the velocity field can be assumed 
divergence-free, divii' = 0. Note that even a slow motion in a 
stratified fluid can be divergent if its spatial scale in radial di- 
rection is not small compared to scale height. In the short-wave 
approximation (in radius) the divergency can, however, be ne- 
glected. 

Our next assumption concerns the pressure. More precisely, 
local thermal disturbances occur at constant pressure so that 
p' lp = -T'/T or s' - -C p p' I p. This assumption is again justi- 
fied by the incompressible nature of the perturbations. Another 
interpretation of this assumption is given by Acheson dl978l l. 
Acheson assumed zero disturbances of total (including mag- 
netic) pressure to involve magnetic buoyancy instability. In our 
derivations, assumptions on constant total or only gas pressure 
are identical because the effect of magnetic buoyancy appears in 
higher order in {krY 1 compared to the terms kept in the equa- 
tions of the next section. 



2.2. Equations 

We start from the linearized equations for small velocity pertur- 
bations, i.e. 



— + (u-V)u' + (u' - V)u + — (V(B-B') 

of /iop 

1 



- (B V) B' - (B' ■ V)B) = - I -VP) + vAh' ( 5 ) 



magnetic field, 
8B' 

= V x (u x B' + u' x B - ifJ x B') . 

8t 



(6) 
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and entropy, 

— +u ■ Vs' +u' ■ Vs 

at 



AT'. 



(7) 



The basic flow is a rotation with nonuniform angular velocity, 
Q = Q(r, 6), and the mean magnetic field is the toroidal one of 
the form (Q]). Perturbations of velocity and magnetic field are ex- 
pressed in terms of scalar potentials after (fJJ and (O. The iden- 
tities 

r (r • V X B') = lT m , r (r • B) = £P m , 

r (r ■ V x u') = iX v , ? (r ■ V x V x u') = it + r 2 ^-)tP v 

are used to reformulate the equations in terms of the potentials. 
The radial component of Eq. ^ then gives the equation for 
the poloidal magnetic field and the radial components of curled 
equations (0 and (0 give the equations for toroidal flow and 
toroidal magnetic field, respectively 

The perturbations are considered as Fourier modes in time, 
azimuth and radius in the form of exp (\(-cjt + m<p + kr)). For an 
instability, the eigenvalue a> should possess a positive imaginary 
part. Only the highest-order terms in kr for the same variable 
were kept in the same equation. 

When deriving the poloidal flow equation, the pressure term 
was transformed as follows 

r-VxVx|-Vpj = -r • V x |-^(Vp) x (VP)| = 

■ V x l-(Vs) x (VP)) = ~ . v X (g x W) = Jr&- 
C p \p ) Cp rC p 

In order to use normalized variables the time is measured in 
units of (Qo is a characteristic angular velocity), the veloc- 
ities are scaled in units of rQo, the normalized frequency (<£>) is 
measured in units of Qo, and the other normalized variables are 
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Introducing the factors kr in the normalizations of poloidal po- 
tentials makes them equal by order of magnitude to the toroidal 
potentials, while in Eqs. © and (|4]i it was P v <k rT v , P m <& rT m . 
The equation for the poloidal flow then reads 



*(iv) = -V(is)-f 2 (iv) 
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1 the radial component of Eq. (f5J curled twice provides the equation 
for the poloidal flow 



where p — cos 8, and 
1= N 



Clftkr 



(10) 



can be understood as special normalization for radial wave- 
length. The first term on the righthand side describes the stabiliz- 
ing effect of the stratification. It vanishes for small A. Apart from 
this stabilizing buoyancy term, the wavelength is only present in 
diffusive terms. The second term of the RHS includes the action 
of finite viscosity, 

vN 2 

e ' = nF' "» 

Similarly, we use below 
_ ?]N 2 _ xN 2 

' ~ e * ~ W 2 



qV 



(12) 



for the diffusive parameters 77 and x- The second and the follow- 
ing lines of Eq. (0 describe the influences of the basic rotation 
and the toroidal field. Only latitudinal derivatives of Q and Qa 
appear. All radial derivatives are absorbed by disturbances which 
vary on much shorter radial scales than Q or Q A . The complete 
system of five equations also includes the equation for toroidal 
flow 

a)(tw) = -i^(tw) + m£l(ZW) - mCl A (tB) 



dp 2 
d_ 

dpW *~ ) ) V~~* ) Qjj 



the equation for toroidal magnetic field 
Cj(£b) = -i^(XB) + wX(qb) - mt(Cl A W) 



cii^ d 1/ 2\ z °" un 
dp y ' dp dp 



+ m 



dp 
,dCl A 
dp 



X 2 dQ. dA 



2 dQ A dV 
dp dp)' 



(14) 



the poloidal field equation 

cj(£a) = -i^(XA)+mQ(XA)-mQ A (XF) (15) 
and the entropy equation, 

ojS = -i^S + milS + LV. (16) 

In the simplest case of Q = Qa = and vanishing diffusivi- 
ties the only nontrivial solution of the above equations are short 
(in radius) gravity waves with 



N 



= — J 1(1 + 1), 1= 1,2,... 

kr 



(17) 



The instabilities we shall find are thus due to either magnetic 
field or nonuniform rotation. It should be kept in mind that the 
above equations are only valid for kr » 1. We shall see that 
the maximum growth rates do indeed belong to the short radial 
scales. 
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2.3. 2D approximation (A » I) 

Generally, the ratio of N/CIq in radiative zones is so large (Fig.Q~|i 
that A can also be large in spite of kr » 1 . Equation (0 in the 
leading order in A then gives S - 0. Then Eqs. ( [To*] ) and ( TT~5T > 
successively yield V — and A = 0. Diffusive terms can also 
be neglected. The equation system reduces to two coupled equa- 
tions for toroidal magnetic field and toroidal flow (Gilman & Fox 
[19971 . 
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(£w) = mn{£W)-mCl A (£B) 

- m ^(( 1 -^) A ) +mB ^(( 1 -^)^)' 
u(Zb) = m£(£lB) - m£(£l A w) . 



(18) 



Here the wave number k drops out and the equations describe 
2D disturbances within decoupled spherical shells. 

The particular case where Q and Qa are both constant sim- 
plifies equations ([T8l strongly. The equations have constant co- 
efficients in this case and can be solved analytically. The solu- 



tion in terms of Legendre polynomials W, B 
eigenfrequencies 



Pf Qi) leads to the 



cl) 1 
m ~ 1(1 + 1) 




1 



1(1+1) J l 2 (l+l) 2 



(19) 



of magnetically modified r-modes (Longuet-Higgins 119641 
Papaloizou & Pringle [T978l ) describing stable horizontal patterns 
drifting in longitude. We shall see in Sec fion lXTl that the case of 
constant Q and Qa shows Tayler instability with 3D approach. 
2D approximation, however, misses the instability. Though the 
result is obtained for particular case of constant Qa it is most 
probably valid in general. 2D incomressive distortions on spher- 
ical surfaces do not change area encircled by toroidal magnetic 
field lines. The circular lines of background field have mini- 
mum length for given encircled area. Any distortion increases 
the length of closed field lines thus increasing magnetic energy. 
There is no possibility to feed an instability by magnetic energy 
release. Only with differential rotation 2D instabilities are possi- 
ble (Gilman & Fox [T997l . 

3. Results and discussion 

We proceed by discussing numerical solutions of the perturba- 
tion equations for special profiles of Qa- In the present paper 
the interaction of toroidal magnetic fields and the basic rotation 
is formulated only for rigid rotation. The work with differential 
rotation is much more complicated, but it is in progress. In the 
Appendix as a first announcement the hydrodynamic stability of 
latitudinal differential rotation in 3D is discussed. 

For the toroidal field two simple geometries are considered. 
First, the quantity Qa is taken constant so that the toroidal field 
has only one belt symmetric with respect to the equator. Second, 
two magnetic belts are considered with equatorial antisymmetry, 
i.e. with a node of B$ at the equator. 

3.1. Fields with equatorial symmetry 

Constant Q and Qa give the simplest realization of the Tayler 
instability. The instability criterion for nonaxisymmetric distur- 
bances (Goossens et al- 1 198 It reads 



2fi 2 - m z 
MM' 



It is satisfied with constant Qa for m = 1 and close to the poles. 

From Section 2.3 we know that the instability can appear 
only in 3D formulation when radial displacements are allowed. 
For this case the equations ©, (fT3l-([l"6l have been solved nu- 
merically. 




2 

N/(Qkr) 

Fig. 2. The stability map for constant Qa and zero diffusivities. 
There is no instability for weak fields with Qa < Q. The thresh- 
old field strength for the instability increases with increasing ver- 
tical wavelength A. 



Consider first ideal fluids with zero diffusivities,^ = rj — v — 
0. Figure [2]gives the resulting stability map. Standard notations, 
Sm and Am, are used for the symmetry types of the global modes 
with azimuthal wave number m symmetric and antisymmetric to 
the equator. However, the disturbances of different fields of the 
same mode possess different symmetries. Rather arbitrarily, our 
symmetry notations correspond to the symmetry of the potential 
W of toroidal field. We did not find instabilities for m + 1 . 

The smallest field strength producing instability corresponds 
to shortest radial scales, X — > 0. The radial velocity is also zero 
in this limit. Note that putting u' r = regardless of which value 
the vertical wavelength, A, has returns to the 2D case without any 
instability. The point here is that for u' r — the entire poloidal 
flow vanishes. The horizontal part of the flow is proportional to 
u' r r/A. A horizontal poloidal flow is necessary for the pinch-type 
instability and the flow remains finite when radial velocity and 
radial scale both approach zero keeping their ratio constant. The 
largest growth rates correspond to this limit. In this sense, the in- 
stability indeed proceeds via 'almost horizontal' displacements. 
The tendency for the instability to prefer indefinitely short radial 
scales shows that finite diffusivities should be included. 

Figure [2] does not show any instability for weak fields with 
Qa < Q. With other words, the basic rotation makes a strong 
stabilizing effect. This may be a special property of our model 
where the ratio Qa/Q is uniform (Pitts & Tayler [19851 I. We in- 
deed find that the threshold field for Qa ~ cos 9 is about ten 
times smaller compared to constant Qa- A much stronger desta- 
bilization, however, is produced by finite diffusion. 

From now on the values 



for n > 0. 
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which are characteristic for the upper part of the solar radiative 
core, are used for the diffusion parameters (fTTb . The relations 
X » J] » v of Eq. d2"TT i are quite typical of stellar radiation 
zones. 
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Fig. 3. The same as in Fig.[2]but with the finite diffusivities (fSTJ. 
The critical field strengths for the onset of the instability are 
strongly reduced compared to ideal MHD. 



For ideal fluids the Tayler instability operates with extremely 
small radial scales. When there is no stabilizing stratification, 
however, finite vertical scales are preferred (Arlt et al. 2007). 
The thermal conductivity decreases the stabilizing effect of the 
stratification and reduces strongly the critical field strengths for 
the instability. The characteristic minima in Fig.[3]correspond to 
small, 1, but finite vertical scales. 

The growth rates for weak fields (Qa < Q) where the in- 
stability exists due to finite diffusion, are very small. For strong 
fields (Qa > Q) the basic rotation is not important and Qa scales 
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Fig. 4. The normalized growth rate <x as function of the magnetic 
field amplitude for the S 1 mode and for A = 0.6 (where the line 
°f Fig EJhas a minimum). The dotted lines for strong, Qa > Q, 
and weak, Q A < Q fields give approximations by the power 
laws o~ ~ Qa and <x ~ Cl 2 A respectively. 
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the growth rate. The dependence of Fig.|4]does indeed approach 
the cr ~ Q A relation (<x is the growth rate) in the strong-field 
limit. The growth rate drops by almost four orders of magni- 
tude when Qa is reduced below 1 . In the weak-field regime it is 
a- ~ Q^ /Q (Spruit fT999l l. 




100 200 300 

Longitude 



Fig. 5. Toroidal field lines of the most rapidly growing eigen- 
modes of Fig.|4]for weak background toroidal field, Qa = 0.1Q 
(upper panel), and strong field, Qa = 10Q (lower panel). The 
horizontal axis gives the longitude. 

The structures of the unstable modes differ between strong 
and weak field regimes. For the weak fields the resulting pattern 
is distributed over the entire sphere. For strong fields it is much 
more concentrated to the poles (Fig [5]) but it remains global in 
latitude. 

3.2. Fields with equatorial antisymmetry 

Now a simple toroidal field with two belts of opposite polarity, 
i.e. 

Q A = bcos6», (22) 

is considered with the diffusivity set (fJTJ. Such a field geome- 
try can result from the action of differential rotation on dipolar 
poloidal fields. Figure [6] shows the corresponding stability map. 
It is similar to the map of Fig. [3] Instability is again found only 
for m — 1 . 

In physical values for the upper part of solar radiation zone 
(with p =i 0.2 g/cm 3 ) one finds 

B a \0 5 b G (23) 
and 

A a 10A Mm, (24) 

so that from Fig.|6]follows a critical magnetic field for the insta- 
bility slightly below 600 G. The mode which first becomes un- 
stable if the field exceeds this critical value has a vertical wave- 
length between 1 and 2 Mm. For larger field strengths, of course, 
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Fig. 6. Stability map for the field model d22l) with two belts and 
equatorial antisymmetry. 
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Fig. 8. The growth rate as function of the toroidal field amplitude 
for A = 0.1. The dotted line shows the parabolic approximation 
& = 0.1b 2 . The scale on the right gives the estimated radial dif- 
fusivity of chemical species after d25l l. 



there is a range of unstable wavelengths. The maximum growth 
rates appears, however, at wavelengths S 1 Mm (Fig.|7j. 




N/{Qkr) N/{Qhr) 



Fig. 7. Growth rates in units of CI for b = 0.01 (left) and b - OA 
(right). The largest rates exist for ^ S 0.1 independent of the 
magnetic field amplitude. All for SI modes. 

The flow field of any instability mixes chemical species also 
in radial direction. Such an instability can thus be relevant to the 
radial transport of the light elements (Barnes et al. 1999). The 
effective diffusivity, Dj u't (u' and I are rms velocity and 
correlation length in radial direction) can roughly be estimated 
from our linear computations assuming that <x m Iju' and i =s 
A/2. With Eq.(E§ this yields 

D T -7- 10 9 & cmV 1 , (25) 

where & is the normalized growth rate given in Fig. [8] For the 
range of 0.01 < b < 0.2 where the plot is closely approximated 
by the parabolic law <x 0.1 b 2 , Eq. (|25l > can be rewritten in 
terms of (Eq. d23"T >) as 

Ct-7.10 4 ^) 2 cmV 1 . (26) 

As known, diffusivities exceeding 10 3 cm 2 s~ 1 in the upper 
radiative core are not compatible with the observed solar lithium 
abundance. Hence, the toroidal field amplitude can only slightly 
exceed the marginal value of about 600 G. In our (simplified) 
formulation, the observed solar lithium abundance seem to ex- 
clud^ any concept of hydromagnetic dynamos driven by Tayler 

2 or rise new problems 



instability in the upper radiation zone of the Sun. We should not 
forget, however, that the superrotation (d£l/dr > 0) at the bottom 
of the convection zone in the equatorial region acts stabilizing 
so that the critical field amplitutes for Tayler instability may be 
higher than the computed 600 G. 

4. Summary 

Linear stability of toroidal magnetic field in rotating stellar ra- 
diation zones is analyzed assuming that the vertical scale of 
the fluctuations is short compared to the local radius ('short- 
wave approximation'). The analysis is global in horizontal di- 
mensions. Stability computations confirm that the most rapidly 
growing perturbations have short radial scales: kr ~ 10 3 . 

We have shown that pinch-type instability of toroidal field 
require nonvanishing radial displacements. The instability does 
not appear in 2D approximation with zero radial velocities. The 
maximum amplitude of stable toroidal magnetic fields for the 
Sun which we have found is about 600 G. This value results only 
for rigid rotation. It will most probably increase if the stabilizing 
influence of the positive radial gradient of Q. in the equatorial 
region of the tachocline is included into the model. 

The field strength in the upper part of the solar radiative 
interior can only marginally exceed the resulting critical val- 
ues. Otherwise the instability would produce too strong radial 
mixing of light elements. After our results all the axisymmet- 
ric hydromagnetic models of the solar tachocline (Riidiger & 
Kitchati nov[T997l Garaud |2002l Sule et al. |2005t Kitchatinov & 
Riidiger 2006) have stable toroidal fields. On the contrary, the 
strong fields ~ 10 5 G which are able to modify g-modes or even 
stronger fields which may influence solar neutrinos are strongly 
unstable with e-folding times shorter than one rotation period. 

3D computations of joint instabilities of toroidal fields and 
differential rotation (Gilman & Fox [19971 Cally |2003l Riidiger et 
al. 2007) can be a perspective for further work. Another tempting 
extension is the inclusion of the poloidal field. The field can be 
important in view of the very short vertical scales of the unstable 
modes. 

In the Appendix we present a calculation with the same equa- 
tions for the hydrodynamic instability of latitududinal differen- 
tial rotation. This instability can already be found in 2D approx- 
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imations (Watson 1981) but it is substantially modified in the 3D 
theory. 
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We see that short rather than long radial scales are preferred. The 
minimum a appears for A ^ 0.6, so that the characteristic wave- 
length of A ^ 6 Mm results after Eq. d24b for the solar tachocline. 
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Fig. A.l. Stability map for the hydrodynamic instability of lati- 
tudinal differential rotation. Most unstable are the perturbations 
of Al symmetry type with the vertical scale A ^ 0.6. The critical 
magnitude of latitudinal shear is reduced to 0.2 1 compared to the 
0.28 value of 2D theory. 



Appendix A: The Watson problem in 3D 

Latitudinal differential rotation can be unstable even without 
magnetic field if the shear 8Q./d8 is sufficiently large (Watson 
1 198 1 1; Dziembowski & Kosovichev 1 1987I I. The instability may 
reduce the differential rotation to its critical value (Garaud 2001) 
which can be relevant to the theory of the solar tachocline. The 
critical relative value of 28% for differential rotation found by 
Watson resulted from a 2D theory (cf. Section [33). The value 
has also appeared in a 3D numerical probe of marginal stability 
of a shell rotating fast enough with the rotation law 

Q = Qo(l -flcos 2 ^), (A.l) 

but for not stratified material (Arlt, Sule & Riidiger 2007). The 
critical shear increases to much higher values, however, if the 
real rotation law (including its radial variation) of the solar 
tachocline is adopted. 

Equations (|9), (Q3), and ([16) of the Section 13721 (in their hy- 
drodynamical version) can also be applied to extend the Watson 
approach by allowance for radial displacements. Note that the 
Reynolds number Re = Q.r 2 /v can be written as 

N 2 /Q. 2 

Re = — - — , (A.2) 

which with Eq. d3T) gives a very large value, O(10 15 ). 

For positive and sufficiently large a, the modes Al and S2 
become unstable. Figure [AT] shows the dependence of the crit- 
ical values of a on the normalized wavelength A (QJJ|. For large 
enough radial wavelengths the 28%-value of the 2D theory is re- 
produced. It is reduced, however, to a = 0.21 in 3D calculations. 



